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Abstract 

We present the preliminary measurement of CP-violating asymmetries in — > {pir)^ tt~^tt~tt^ 
decays using a time-dependent Dalitz plot analysis. The results are obtained from a data sample 
of 213 million T{4S) — > BB decays, collected by the BABAR detector at the PEP-II asymmetric- 
energy B Factory at SLAC. This analysis extends the narrow-p quasi-two-body approximation 
used in the previous analysis, by taking into account the interference between the p resonances of 
the three charges. We measure 16 coefficients of the bilinear form factor terms occurring in the 
time-dependent decay rate of the B^ meson with the use of a maximum-likelihood fit. We derive 
the physically relevant quantities from these coefficients. We measure the direct CP-violation 
parameters Ap^ = -0.088 ± 0.049 ± 0.013 and C = 0.34 ± 0.11 ± 0.05, where the first errors are 
statistical and the second systematic. For the mixing-induced CP-violation parameter we find 
S = —0.10 lb 0.14 lb 0.04, and for the dilution and strong phase shift parameters respectively, we 
obtain AC = 0.15 ib 0.11 ib 0.03 and AS = 0.22 ib 0.15 ib 0.03. For the angle a of the Unitarity 
Triangle we measure (1131^7 =t 6)°, while only a weak constraint is achieved at the significance 

level of more than two standard deviations. Finally, for the relative strong phase 5^ between the 

B^ — > p~7r~^ and B^ — > p+Tr" transitions we find (— 67^3^ ib 7)°, with a similarly weak constraint 
at two standard deviations and beyond. 
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1 INTRODUCTION 



Measurements of the parameter sin2/3 [1,2] have estabhshed CP violation in the meson system 
and provide strong support for the Kobayashi and Maskawa model of this phenomenon as arising 
from a single phase in the three-generation CKM quark-mixing matrix [3] . We present in this letter 
preliminary results from a time-dependent analysis of the — 7r+7r~7r° Dalitz plot (DP) that is 
dominated by the p intermediate resonances. The goal of the analysis is the simultaneous extraction 
of the strong transition amplitudes and the weak interaction phase a = arg [— ^(i^^/V^^V^,] of the 
Unitarity Triangle. In the Standard Model, a non-zero value for a would be responsible for the 
occurrence of mixing-induced CP violation in this decay. The BABMi and Belle experiments have 
obtained constraints on a from the measurement of effective quantities sin2Q;eff in B decays to 
vr+vr^ [4,5] and BABAR from B decays to [6], using an isospin analysis [7] that involves the 

other charges of these final states to obtain bounds on a — aeff • 

Unlike 7r"^7r~, p^Tr^ is not a CP eigenstate, and four flavor-charge configurations {B^{B^) — >■ 
p^TT^) must be considered. The corresponding isospin analysis [8] is unfruitful with the present 
statistics since two pentagonal amplitude relations with 12 unknowns have to be solved (compared 
to 6 unknowns for the tt+tt" and p'^ p~ systems). However, it has been pointed out by Snyder and 
Quinn [9] , that one can obtain the necessary degrees of freedom to constrain a without ambiguity by 
explicitly including in the analysis the variation of the strong phases of the interfering p resonances 
in the Dalitz plot. 

The present analysis focuses on the decay B^ — > (/ott)" — > 7r+7r~7r°. The data sample used is a 
superset of that used for our previous time-dependent result [10] , obtained within the quasi- two- 
body approximation. In that approach, the analysis was restricted to the chargcd-p regions in the 
Dalitz plot, and the interference regions were removed. Here, we extend the analysis to the entire 
region of interest in the Dalitz plot, which contains the p resonances of all three charges and their 
interference. 



1.1 DECAY AMPLITUDES 

We consider the decay of a spin-zero B^ with four-momentum ps into the three daughters 7r^(p+), 
^ (P-)) ^^(Po)) with corresponding four- momenta. Using as independent (Mandelstam) variables 
the invariant masses-squared 

s+ = {p++Paf , s- = {p-+Pof , (1) 

the invariant mass of the positive and negative pion, sq = {p+ -|-p_)^, is obtained from energy 
and momentum conservation 

So = tn^BO + 2m^+ + m^o - s+ - S- . (2) 

The differential B^ decay width with respect to the variables defined in Eq. (1) {i.e., the Dalitz 
plot) reads 

dr(B«^.W)= pL3^rf»,d,_, (3) 

where A3Tr is the Lorentz-invariant amplitude of the three-body decay. 

We assume in the following that the amplitudes ^3,^ and its complex conjugate A^t^, corre- 
sponding to the transitions B^ — > 7r+7r~7r° and B^ 7r"'"7r~7r'', respectively, are dominated by the 
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three resonances p"*", p~ and p . The p resonances are assumed to be the sum of the ground state 
p(770) and the radial excitations p{lA5Q) and p(1700), with an initial set of resonance parame- 
ters and relative amplitudes determined by a combined fit to t'^ — > Ft-tt^tt^ and e^e" — tt'^tt" 
data [11]. Since the hadronic environment is different in B decays, we cannot rely on this result 
and therefore determine the relative /3(1450) amplitude simultaneously with the CP parameters 
from the fit. Variations of the other parameters and possible contributions to the 7r"'"7r~7r° 
decay other than the p's are studied as part of the systematic uncertainties (Section 4). 

Including the B^B^ mixing parameter q/p into the B^ decay amplitudes, we can write [9,12] 

As^ = f+A+ + f_A- + foA\ (4) 
As^ = UA+ + f_A- + foA\ (5) 

where the (with the p charge k = {+, — , 0}) are functions of the Dalitz variables s+ and s_ that 
incorporate the kinematic and dynamical properties of the B^ decay into a (vector) p resonance 
and a (pscudoscalar) pion, and where the A'^ are complex amplitudes^ that may comprise weak 
and strong transition phases and that arc independent of the Dalitz variables. Note that the 
definitions (4) and (5) imply the assumption that the relative phases between the p(770) and its 
radial excitations are CP-conserving. 

Following Rcf. [11], the p resonances arc parameterized by a modified relativistic Breit-Wigner 
function introduced by Gounaris and Sakurai (GS) [15]. Due to angular momentum conservation, 
the spin-one p resonance is polarized in a helicity-zero state. For a p'^ resonance with charge k, 
the GS function is multiplied by the kinematic function — 4|p^||pt-| cos^^j where the momenta are 
defined in the p-resonancc rest frame, and where Pr is the momentum of the particle not from p 
decay, and cos9^ the cosine of the helicity angle of the p'^. For the p~^ (p ), 0+ (^-) is defined by 
the angle between the vr" (tt") in the p"*" (p~) rest frame and the p"*" (p~) flight direction in the B^ 
rest frame. For the p", Oq is defined by the angle between the tt"^ in the p° rest frame and the p° 
flight direction in the B'^ rest frame. With these definitions, each pair of GS functions interferes 
destructively at equal masses-squared. 

The occurrence of cos in the kinematic functions substantially enhances the interference 
between the different p bands in the Dalitz plot, and thus increases the sensitivity of this analysis [9]. 



1.2 TIME DEPENDENCE 



With At = tsn — hag defined as the proper time interval between the decay of the fully reconstructed 
B^^ and that of the other meson B^^^, the time-dependent decay rate |.4.3^(At)|^ (|^3j^(At)|^) when 
the tagging meson is a B^ (-B") is given by 



g-IAtl/r^o 



=F 

± 2Im 



^SttI^ - I^SttI^) cos(AmdAi) 



AznAl^ sm{AmdAt) 



(6) 



where r^o is the mean B^ lifetime, Arrid is the B^B^ oscillation frequency, and where we have 
assumed that CP violation in B^B^ mixing is absent {\q/p\ = 1), AF^^ = and CPT is conserved. 

®The amplitude superscript "/t" denotes the charge of the p from the decay of the B° meson. 
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Inserting the amplitudes (4) and (5), one obtains for the terms in Eq. (6) 



|^3,r|'±|^37r|' 
Im (^Stt^Stt) 



Y: \U\'ut + 2 Y: ( Re [/«/*] C/«1;^^-Im [/«/*] , 

k;6{+,-,0} K<(7e{+,-,0} 

Ke{+,-,o} K<ae{+,-,o} 



with 



jj±,Re(Ini) 



rRe 

KIT 

7-Im 



Re(Im) [A'^A'^* ± l"!'"* 
Im [A'^A''* 
Re [A'^A''* - A" A'' 
Im [I'^A'^* + A" A" 



(8) 
(9) 
(10) 

(11) 
(12) 



The 27 coefficients (8)-(12) are real- valued parameters that multiply the fnf^ bilinears [16]. 
They are the observables that are determined by the fit. Each of the coefficients is related in a 
unique way to the physically more intuitive quantities, like tree-level and penguin-type amplitudes, 
the angle a, or the quasi-two-body CP and dilution parameters [10] (c/. Section 6). 

The coefficients are related to resonance fractions (branching fractions and charge asymme- 
tries), the U~ determine the relative abundance of the decay into p'^7r~ and p^ir^ (dilution) 
and the time-dependent direct CP asymmetries. The 1^ measure mixing-induced CP violation and 
are sensitive to strong phase shifts. Finally, the uta^'^^^"^^ and I^a^^"^^ describe the interference 
pattern in the Dalitz plot. Their presence distinguishes this analysis from the previous quasi-two- 
body one [10]. They represent the additional degrees of freedom that allows one to determine the 
unknown penguin pollution and the relative strong phases. However, because the overlap regions 
of the resonances are small and because the event reconstruction in these regions suffers from large 
misreconstruction rates and background, a substantial data sample is needed to perform a fit that 
constrains all amplitude parameters. 

The choice to fit for the U and / coefficients rather than fitting for the complex transition 
amplitudes and the weak phase a directly is motivated by the following technical simplifications: 
(i) in contrast to the amplitudes, there is a unique solution for the U and / coefficients requiring only 
a single fit to the selected data sample^, (ii) in the presence of background, the U and / coefficients 
are approximately Gaussian distributed, which in general is not the case for the amplitudes, and 
{in) the propagation of systematic uncertainties and the averaging between different measurements 
are straightforward for the t/'s and /'s. 

We determine the quantities of interest in a subsequent least-squares fit to the measured U and 
/ coefficients. 



'^The parameterization (7) is general: the information on the mirror solutions {e.g., on the angle a) that are 
present in the transition amplitudes A'^, A"' is conserved. 
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1.3 NORMALIZATION 



The decay rate (6) is used as probability density function (PDF) in a maximum-likelihood fit and 
must therefore be normalized: 

/M 12! IT 12x 1-^3^(^^)1' , (13) 

where 

{\As.\' + \As.\') = E + 2Re E (A/;) (^.^/^ + ^ • ^.Vi • (14) 

Ke{+,-,Q} K<ae{+,-,o} 

The complex expectation values {f^fa) obtained from high-statistics Monte Carlo integration 
of the Dalitz plot (3), taking into account acceptance and resolution effects. 

The normalization of the decay rate (6) renders the normalization of the U and / coefficients 
arbitrary, so that we can fix one coefficient. By convention, we set = 1. 

For a small — ^ p^i^^ signal, the time-dependent CP information from this mode is marginal. 
As a consequence, we are allowed to simplify the model by setting the coefficients of the sine and 
cosine terms that involve this mode to zero, which entails only a small systematic error on the 
result. This reduces the number of free parameters from 26 to 16. 

1.4 THE SQUARE DALITZ PLOT 

The signal events and also the combinatorial e^e" qq {q = u, d, s, c) continuum background 
events populate the kinematic boundaries of the Dalitz plot due to the low final state masses 
compared to the mass. Hence the representation Eq. (3) is inadequate when one wants to 
use empirical reference shapes in a maximum- likelihood fit. Another practical disadvantage of the 
Dalitz variables is that they are not constant along the kinematic boundaries of the B^ — tt+tt^tt'^ 
Dalitz plot, which complicates the correction of efficiency and resolution effects. We therefore apply 
the transformation 

ds^ ds^ — > I dct J| dm dO' , (15) 
which defines the Square Dalitz plot (SDP). The new coordinates are 

m' = - arccos { 2 "^0 " "^0"" _ 1^ and 6' = \ , (16) 

where mo is the invariant mass between the charged tracks, m^'^^ = m^o — m^o and m™™ = 2m^+ 
are the kinematic limits of mo, 9o is the helicity angle, and J is the Jacobian of the transformation 
that zooms into the kinematic boundaries of the Dalitz plot. The new variables have validity ranges 
between and 1. The determinant of the Jacobian is given by 

1^ + 71 ^ I * II *i ^"^0 5cos6lo , . 

|detJ| = 4|p+llpo|mo--^--^ , (17) 

where |p+| = ^ E'^ — and |Pol = \J ~ '^^o 1 and where the energies E*^. and E^ are in the 
7r+7r~ rest frame. 

Figure 1 shows the original (left hand plot) and the transformed (right hand plot) Dalitz plots for 
Monte Carlo B^ tt+tt^tt^ events generated according to Eqs. (4) and (5) with equal abundance 
of all p charges and with vanishing relative strong phase. The plots illustrate the homogenization 
of the Dalitz plot obtained after the transformation (15). 
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Figure 1: Nominal (left) and square (right) ■k~^7:~t:^ Dalitz plots obtained from Monte Carlo 

generated events without detector simulation. The generating amplitudes used are A'^ = A~ = 
= 1 so that they interfere destructively at equal p masses. The hatched areas indicate the 
main overlap regions between the different p bands. The contour lines in both plots correspond to 
y^sIjI~~o = 1.5 GeV/c^. The middle plot depicts the Jacobian determinant (17) of the transforma- 
tion (15). The plot shows the distribution one would obtain in the square Dalitz plot for uniformly 
distributed events in the nominal Dalitz plot. 

2 THE BABAR DETECTOR AND DATASET 

The data used in this analysis were collected with the BABAR detector at the PEP-II asymmetric- 
energy e~^e~ storage ring at SLAC between October 1999 and July 2004. The sample consists of 
about 192 fb~\ corresponding to (213.1±2.3) x 10^ BB pairs collected at the T(45) resonance ("on- 
resonance"), and an integrated luminosity of 11.6 fb~^ collected about 40 MeV below the T{4S) 
( "off-resonance" ) . 

A detailed description of the BABAR detector is presented in Ref. [17]. The tracking system used 
for track and vertex reconstruction has two main components: a silicon vertex tracker (SVT) and a 
drift chamber (DCH), both operating within a 1.5 T magnetic field generated by a superconducting 
solenoidal magnet. Photons are identified in an electromagnetic calorimeter (EMC) surrounding 
a detector of internally reflected Cherenkov light (DIRC), which associates Cherenkov photons 
with tracks for particle identification (PID). Muon candidates are identified with the use of the 
instrumented flux return (IFR) of the solenoid. 

3 ANALYSIS METHOD 

The U and / coefficients and the B^ tt~^'k~tt^ event yield are determined by a maximum- 
likelihood fit of the signal model to the selected candidate events. Kinematic and event shape 
variables exploiting the characteristic properties of the events are used in the fit to discriminate 
signal from background. We limit the size of the data sample that enters the fit by tightening the 
acceptance requirements for the discriminant variables compared to similar analyses, because the 
fit model with at least 17 physical parameters is rather involved. With the same goal and because 
the modeling of the distribution of the continuum events in the Dalitz plot is delicate, we remove 
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the center of the Dahtz plot from the analysis. This requirement does not affect the signal. 

3.1 EVENT SELECTION AND BACKGROUND SUPPRESSION 

We reconstruct — >■ tt+tt^tt^ candidates from pairs of oppositely-charged tracks, forming a good 
quality vertex, and a 7r° candidate. We use information from the tracking system, EMC, and DIRC 
to remove tracks for which the PID is consistent with the electron, kaon, or proton hypotheses. In 
addition, we require that at least one track has a signature in the IFR that is inconsistent with the 
muon hypothesis. The 7r° candidate mass must satisfy 0.11 < 771(77) < 0.16GeV/c^, where each 
photon is required to have an energy greater than 50 MeV in the laboratory frame (LAB) and to 
exhibit a lateral profile of energy deposition in the EMC consistent with an electromagnetic shower. 

Two requirements are applied on the Dalitz plot. Firstly, the invariant mass of the two tracks mo 
must be larger than 0.52GeV/c2. This rejects about 80% of the B° ^ K^{-^ TT'^TT )7r background 
events, which due to the long lifetime of the would require a dedicated At treatment. This cut 
retains 98% (100%) of signal — > p^TX^ {B^ p^vr^) events. Secondly, we remove the center of 
the Dalitz plot by requiring that at least one of the three invariant masses, mo, m+ or m_, is lower 
than 1.5 GeV/c^. 

A S-meson candidate is characterized kinematically by the energy-substituted mass mss = 
[{^s + po • ps)^/£'o — p^]2 and energy difference = E*^ — \^/s, where {Eb,Pb) and {Eo,pq) 
are the four-vectors of the S-candidate and the initial electron-positron system, respectively. The 
asterisk denotes the T(4S) frame, and s is the square of the invariant mass of the electron-positron 
system. We require 5.272 < rriES < 5.288 GeV/c^, which retains 81% of the signal and 8% of the 
continuum background events. The AE resolution exhibits a dependence on the vr'^ energy and 
therefore varies across the Dalitz plot. Wc account for this effect by introducing the transformed 
quantity A^;' = {2AE-AE+-AE_)/{AE+-AE_), with AE±{mo) = c±-(c± ^ c) (mo/mS'^^)^, 
where mo monitors the 7r°-energy dependence. We use the values c = 0.045 GeV, c_ = —0.140 GeV, 
c+ = 0.080 GeV, mg"^^ = 5.0 GeV, and require —1 < AE' < 1. These settings have been ob- 
tained from Monte Carlo (MC) simulation and are tuned to maximize the selection of correctly 
reconstructed over misreconstructed signal events. The cut retains 75% (25%) of the signal (con- 
tinuum). 

Backgrounds arise primarily from random combinations in continuum events. To enhance dis- 
crimination between signal and continuum, we use a neural network (NN) to combine four discrim- 
inating variables: the angles with respect to the beam axis of the B momentum and B thrust axis 
in the T{4S) frame, and the zeroth and second order monomials Lo,2 of the energy flow about the 
B thrust axis. The monomials are defined by Lj = J2iPi ^ |cos^ip, where 9i is the angle with 
respect to the B thrust axis of track or neutral cluster i, pi is its momentum, and the sum excludes 
the B candidate. The NN is trained in the signal region with off-resonance data and simulated 
signal events. The final sample of signal candidates is selected with a cut on the NN output that 
retains 77% (8%) of the signal (continuum). 

The time difference At is obtained from the measured distance between the z positions (along 
the beam direction) of the B^^ and B^^^ decay vertices, and the boost P'j = 0.56 of the e~^e~ 
system^. To determine the flavor of the B^^^^ we use the tagging algorithm of Ref. [18]. This 
produces four mutually exclusive tagging categories. We also retain untagged events in a fifth 
category to improve the efficiency of the signal selection and because these events contribute to the 
measurement of direct CP violation. Events with multiple B candidates passing the full selection 

^At is defined as: At = Az/f3-/c 
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occur in 16% (p^TT^) and 9% (p^tt^) of the cases. If the multiple candidates have different vr^'s, 
we choose the candidate with the reconstructed 7r° mass closest to the nominal one; if not, one 
candidate is selected at random. 

The signal efficiency determined from MC simulation is 24% for — > p^ir^ and — ^ p'^ir^ 
events, and 11% for non-resonant B^ — >■ 7r"'"7r~7r'^ events. 

Of the selected signal events, 22% (5° P^tt^), 13% (B° y9°7r°), and 6% (non-resonant) 
are misreconstructed, mostly due to combinatorial background from low-momentum tracks and 
photons. They concentrate in the corners of the Dalitz plot. The fraction of misreconstructed 
events strongly varies across the tagging categories. 



3.2 BACKGROUND FROM OTHER B DECAYS 

We use MC-simulated events to study the background from other B decays. The exclusive B- 
background modes are grouped into eighteen classes with similar kinematic and topological prop- 
erties. More than hundred decay channels have been considered of which thirty-six are retained 
in the likelihood model. The most significant ones are B^ P~^P^ with longitudinal polarization 
(27±18 events expected), B+ 7r+p° (48±6), B+ 7r°p+ (43±7), B° p+p' with longitudinal 
polarization (50± 10), B^ (ai7r)° (29± 11), S° ^ p-K+ (61 ± 11), arid S° ^ higher kaon reso- 
nances (6 lb 1). The charmed modes D~(^ 7r~7r°)7r+ and B^ D^{—>^ 7r"'"7r~)7r° contribute 
to the selected data sample and are considered in individual classes. They do not interfere with 
the signal due to the long D lifetime. We also assign classes to the modes B^ D^{-^ K~^7r~)7r^ 
and B'^ — J/il){-^ t^f)'K^. In total we expect 49 it 15 exclusive 6 — c events. Two additional 
classes account for inclusive neutral and charged 6 — c decays, where we expect 82 it 6 and 181 ± 9 
events, respectively. 



3.3 THE MAXIMUM-LIKELIHOOD FIT 

We perform an unbinned extended maximum-likelihood fit to extract the inclusive B^ — > tt^tt^tt^ 
event yield and the U and / coefficients defined in Eqs. (8)-(12). The fit uses the variables tties, 
A£'', the NN output, and the Dalitz plot to discriminate signal from background. The At mea- 
surement allows to determine mixing-induced CP violation and provides additional continuum- 
background rejection. 

The selected on-resonance data sample is assumed to consist of signal, continuum-background 
and 5-background components, separated by the flavor and tagging category of the tag side B 
decay. The signal likelihood consists of the sum of a correctly reconstructed ( "truth- matched" , 
TM) component and a misreconstructed ( "self-cross- feed" , SCF) component. 

The probability density function (PDF) Vf for an event i in tagging category c is the sum of 
the probability densities of all components, namely 



'Pi = ^Snfin (1 - /sCF)'P37r-TM,i + /sCF^Stt-SCF, 
+ ^qq2 + Qtag,iAqq,ta.g) Vqq^i 

class -| 

+ ^ N B+j /b+ J ^5 ( 1 + 9tag,j^B+ , tag J ) ^3+ ,ij 
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+ E ^SO,/^o,no,, , (18) 

where: N^-j^ is the total number of 7r'^7r~7r° signal events the data sample; /^^ is the fraction 
of signal events that are tagged in category c; /g^p is the fraction of SCF events in tagging 
category c, averaged over the Dalitz plot; "P^tt-tm i ^Stt-SCF i products of PDFs of the 

discriminating variables used in tagging category c for TM and SCF events, respectively; N^^ is 
the number of continuum events that are tagged in category c; gtag,i is the tag flavor of the event, 
defined to be +1 for a B^^^ and —1 for a i?tagi ^qq, tug parameterizes possible tag asymmetry in 
continuum events; V^^^ is the continuum PDF for tagging category c; N^^^^ (-^dass) number 
of charged (neutral) i?-related background classes considered in the fit; N^+j (N^Uj) is the number 
of expected events in the charged (neutral) S-background class j; f^+j (/ijcj) the fraction of 
charged (neutral) 5-background events of class j that are tagged in category c; ^s+^tag.j describes 
a possible tag asymmetry in the charged-B background class j; correlations between the tag and 
the position in the Dalitz plot (the "charge") are absorbed in tag-flavor-dependent Dalitz plot 
PDFs that are used for charged-B and continuum background; ^. is the S+-background PDF 
for tagging category c and class j\ flnally, V%,o ■■ is the neutral- S-background PDF for tagging 
category c and class j. 

The PDFs Vx are the product of the four PDFs of the discriminating variables, xi = niEs, 
X2 = AE', xs = NNoutput, and the triplet x^ = {m', 9', At}: 

^AJ) - n^u-)(^^)- (19) 

k=l 

The extended likelihood over all tagging categories is given by 

-c^n^-^'fi^^' (20) 

c=l i 

where A'^'^ is the total number of events expected in category c. 

A total of 39 parameters, including the inclusive signal yield and the parameters from Eq. (6), 
are varied in the fit. 

3.3.1 THE At AND DALITZ PLOT PDFS 

Signal Parameterization. The Dalitz plot PDFs require as input the Dalitz plot-dependent 
relative selection efficiency, e = e{m',9'), and SCF fraction, /scF = fscF{n^',9'). Both quantities 
are taken from MC simulation. They are given in Fig. 2 (left plot for e and right plot for /scf), 
where the symmetry of the Dalitz plot has been used to fold the upper 6' half into the lower one. 
Away from the Dalitz plot corners the efficiency is uniform, while it decreases when approaching the 
corners, where one out of the three bodies in the final state is close to rest so that the acceptance 
requirements on the particle reconstruction become incisive. Combinatorial backgrounds and hence 
SCF fractions are large in the corners of the Dalitz plot due to the presence of soft neutral clusters 
and tracks. 

For an event i, we define the time-dependent Dalitz plot PDFs 

P3^-TM,i = ei{l- f SCF, i)\det Jill Ai^{At)W (21) 

P3n-scF,i = eifscF,i\detJi\\At{At)W (22) 
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Figure 2: Selection efficiency of vr^vr vr*^ events (left) and fraction of misreconstructed events 

(right) in the (symmetrized) square Dalitz plot for MC-simulated events. 



where Pstt-tm,? and -P37r-sCF,i are normalized. The corresponding phase space integration involves 
the expectation values (£(1- fsCF)\det J\ fT*) and (e /scf 1 dct J| /'^/'^*) for TM and SCF 
events, where the indices k, a run over all resonances belonging to the signal model. The expectation 
values are model-dependent and are computed with the use of MC integration over the square Dalitz 
plot: 

_ Jo Jo £ (1 - /scf) I det J| rr* dm'de' 



(e(l-/scF)|detJ| 



(23) 



/oVo£|detJ|/«r*dmW 

and similarly for (e | det J| f^f^*), where all quantities in the integrands are Dalitz plot-dependent. 

Equation (18) invokes the phase space-averaged SCF fraction fscF = (/scf | det J| f^f"*). As 
for the PDF normalization, it is decay-dynamics-dependent. It has to be computed iteratively, 
though the remaining systematic uncertainty after one iteration step is small. We determine the 
average SCF fractions separately for each tagging category from MC simulation. 

The width of the dominant /9(770) resonance is large compared to the mass resolution for TM 
events (about 8MeV/c^ core Gaussian resolution). We can therefore neglect resolution effects in 
the TM model. Misreconstructed events have a poor mass resolution that strongly varies across 
the Dalitz plot. It is described in the fit by a 2 x 2-dimensional resolution function 



(24) 



which represents the probability to reconstruct at the coordinate an event that has the 

true coordinate (mj,^^). It obeys the unitarity condition 



1 1 



j J RscF{m'„ 9'„ m[, 9[) drn'^dO', = 1 , V {m[, e[) G SDP , 



(25) 







and is convolved with the signal model. The i?sCF function is obtained from MC simulation. 
Figure 3 shows the resolution function of TM (left) and SCF events (right) for two coordinates 
depicted by the open stars. 

We use the signal model described in Section 1.1. It contains the dynamical information and 
is connected with via the matrix element (6), which serves as PDF. It is diluted by the effects 
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Figure 3: Resolution for TM (left) and SCF events (right hand plot) in the square Dalitz plot for 
two coordinates indicated by the open stars. 



of mistagging and the limited vertex resolution [10]. The At resolution function for signal and 
S-background events is a sum of three Gaussian distributions, with parameters determined by a 
fit to fully reconstructed decays [18]. 

Background Parameterization. The Dalitz plot- and At-dependent PDFs factorize for the 
charged- -B-background modes, but not (necessarily) for the neutral-i? background due to B^B^ 
mixing. 

The charged 5-background contribution to the likelihood (18) invokes the parameter ^4^+ ^^^g, 
multiplied by the tag flavor gtag of the event. In the presence of significant tag- "charge" correla- 
tion (represented by an effective flavor-tag- versus-Dalitz-coordinate correlation), it parameterizes 
possible direct CP violation in these events. We also use distinct square Dalitz plot PDFs for each 
reconstructed B flavor tag, and a flavor-tag-averaged PDF for untagged events. The PDFs are 
obtained from MC simulation and are described with the use of non-parametric functions. The 
At resolution parameters are determined by a fit to fully reconstructed B'^ decays. For each B^- 
background class we adjust effective lifetimes to account for the misreconstruction of the event that 
modifies the nominal At resolution function. 

The neutral- i? background is parameterized with PDFs that depend on the flavor tag of the 
event. In the case of CP eigenstates, correlations between the flavor tag and the Dalitz coordinate 
are expected to be small. However, non-CP eigenstates, such as a^vr^, may exhibit such correlation. 
Both types of decays can have direct and mixing-induced CP violation. A third type of decays 
involves charged kaons and does not exhibit mixing-induced CP violation, but usually has a strong 
correlation between the flavor tag and the Dalitz plot coordinate (the kaon charge), because it 
consists of i?-flavor eigenstates. The Dalitz plot PDFs are obtained from MC simulation and are 
described with the use of non-parametric functions. For neutral B background, the signal At 
resolution model is assumed. 

The Dalitz plot treatment of the continuum events is similar to the one used for charged- -B 
background. The square Dalitz plot PDF for continuum background is obtained from on-resonance 
events selected in the rriEs sidebands and corrected for feed-through from B decays. A large 
number of cross checks has been performed to ensure the high fidelity of the empirical shape 
parameterization. Analytical models have been found insufficient. The continuum At distribution 
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is parameterized as the sum of three Gaussian distributions with common mean and three distinct 
widths that scale the At per-event error. This yields six shape parameters that are determined by 
the fit. The model is motivated by the observation that the At average is independent of its error, 
and that the At RMS depends linearly on the At error. 

3.3.2 PARAMETERIZATION OF THE OTHER VARIABLES 

The rn-ES distribution of TM signal events is parameterized by a bifurcated Crystal Ball func- 
tion [20], which is a combination of a one-sided Gaussian and a Crystal Ball function. The mean 
of this function is determined by the fit. A non-parametric function is used to describe the SCF 

signal component. 

The AE' distribution of TM events is parameterized by a double Gaussian function, where 
all five parameters depend linearly on uiq. Misreconstructed events are parameterized by a broad 
single Gaussian function. 

Both TTiES and AE' PDFs are parameterized by non-parametric functions for all 5-background 
classes. 

The rriES and AE' PDFs for continuum events are parameterized with an Argus shape func- 
tion [21] and a second order polynomial, respectively, with parameters determined by the fit. 

Wc use non-parametric functions to empirically describe the distributions of the NN outputs 
found in the MC simulation for TM and SCF signal events, and for i?-background events. We 
distinguish tagging categories for TM signal events to account for differences observed in the shapes. 

The continuum NN distribution is parameterized by a third order polynomial that is defined to 
be positive. The coefficients of the polynomial arc determined by the fit. Continuum events exhibit 
a correlation between the Dalitz plot coordinate and the shape of the event that is exploited in the 
NN. The tight requirement that eliminates the center of the Dalitz plot has the purpose to reduce 
such correlation. To correct for residual effects, we introduce a linear dependence of the polynomial 
coefficients on the distance of the Dalitz plot coordinate to the kinematic boundaries of the Dalitz 
plot. The parameters describing this dependence are determined by the fit. 

4 SYSTEMATIC STUDIES 

The contributions to the systematic error on the signal parameters are summarized in Table 1. 

The uncertainties associated with Arrid and r are estimated by varying these parameters within 
the uncertainties on the world average [19]. 

The systematic effects due to the signal PDFs ("Signal description" field in Table 1) comprise 
uncertainties in the PDF parameterization, the treatment of misreconstructed events, the tagging 
performance, and the modeling of the signal contributions. 

When the signal PDFs are determined from fits to a control sample of fully reconstructed B 
decays to exclusive final states with charm, the uncertainties are obtained by varying the parameters 
within the statistical uncertainties. In other cases, the dominant parameters have been left free to 
vary in the fit, and the differences observed in these fits are taken as systematic errors. 

The average fraction of misreconstructed signal events predicted by the MC simulation has been 
verified with fully reconstructed B — Dp events [10]. No significant differences between data and 
the simulation were foTind. Wc vary fscF for all tagging categories by relative 25% to estimate the 
systematic uncertainty (contained in "Signal description" in Table 1). 
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0.034 


0.097 


0.099 


0.019 


0.059 


0.012 


0.028 


Fixing 10 p^ir^ parameters 


0.019 


0.010 


0.024 


0.004 


0.110 


0.069 


0.102 


0.000 


Fit Bias 


0.078 


0.076 


0.221 


0.240 


0.051 


0.048 


0.055 


0.047 


Sum 


0.135 


0.171 


0.704 


0.467 


0.179 


0.149 


0.170 


0.175 



Table 1: Summary of systematic uncertainties. 



Tagging efficiencies, dilutions and biases for signal events are varied within their experimental 
uncertainties. 

The most important contribution to the systematic uncertainty stems from the signal modeling 
of the Dalitz plot dynamics. We vary the mass and width of the p(770) and p(1450) within ranges 
that exceed twice the errors found for these parameters in the fits to r and e+e~ data [11]. Since 
some of the U and / coefficients exhibit significant dependence on the /9(1450) contribution, we leave 
its amplitude (phase and fraction) free to vary in the nominal fit. We vary the relative amount 
of p(1700) by 30% with respect to the nominal model, and its phase by 8°, to assign a systematic 
error. We have performed a fit where the p(1700) amplitude parameters (magnitude and phase) 
are free to vary. We find results that are in agreement with the nominal model. The variations 
for the U and / coefficients observed in this fit compared to the nominal one are smaller than the 
systematic uncertainties we assign due to the p(1700) amplitude uncertainty. 

To estimate the contribution from non-resonance B^ — ^ 7r'^7r~7r° events, we have performed an 
independent analysis where we apply the contrary of the Dalitz plot requirement that is used in 
the nominal analysis; to remove the p signal we retain only those events for which the minimum 
invariant mass exceeds 1.5GeV/c^. For simplicity, we assume a uniform Dalitz distribution for the 
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non-resonance events. The fit finds no non-resonant events so that we can determine a prehminary 
upper Hmit of 1.4 x 10~^ at 90% confidence level (statistical errors only). According to this limit, 
we add simulated non-resonant events to the nominal data sample to estimate the systematic 
uncertainty. We have also searched for the presence of — /o(980)7r° events without finding 
evidence for a signal. 

A major source of systematic uncertainty is the B-background model. The expected event yields 
from the background modes are varied according to the uncertainties in the measured or estimated 
branching fractions. Since S-background modes may exhibit CP violation, the corresponding 
parameters are varied within appropriate uncertainty ranges. As is done for the signal PDFs, we 
vary the At resolution parameters and the flavor-tagging parameters within their uncertainties and 
assign the differences observed in the on-resonance data fit with respect to nominal fit as systematic 
errors. 

The parameters for the continuum events are determined by the fit. No additional systematic 
uncertainties are assigned to them. An exception to this is the Dalitz plot PDF: to estimate 
the systematic uncertainty from the mss sideband extrapolation, we select large samples of off- 
resonance data by loosening the requirements on AE and the NN. We compare the distributions 
of m' and 9' between the tbes sideband and the signal region. No significant differences are found. 
We assign as systematic error the effect seen when weighting the continuum Dalitz plot PDF by 
the ratio of both data sets. This effect is mostly statistical in origin. To account for possible 
inaccuracies in the empirical parameterization, we add test components with floating event yields 
to the fit that consist of continuum-background-like reference distributions in all fit variables but 
the Dalitz plot, for which signal- like distributions are used. In addition, we leave the S-background 
classes free to vary in the fit. This allows the fit to absorb events that, due to possible problems 
in the continuum Dalitz plot description, would bias the signal yield. The shifts in the signal 
parameters observed when adding these test components arc taken as systematic uncertainties 
("continuum parameterization"' in Table 1). It leads to significant effects in the p^n^ region of the 
Dalitz plot. 

We assess the dependence of the results on whether all the 27 U and / coefficients are free in 
the fit or only the 16 most significant ones. To study this effect, we generate MC samples with 
values for the U and / coefficients according to our nominal 16 parameter fit result. We then 
perform a full amplitude fit to these 16 coefficients (c/. Section 6) and compute from the best fit 
the expected values for the missing coefficients. Monte Carlo samples are generated with the use of 
all 27 coefficients, which are fit with the nominal model where 10 coefficients are set to zero. The 
observed systematic effect on the measured U and / coefficients depend on the branching fraction 
for — p^TT^ for which we use our upper bound [22] (which is in agreement with the finding in 
this analysis). 

Finally, to validate the fitting tool, we perform fits on large MC samples with the measured 
proportions of signal, continuum and i?-background events. No significant biases are observed in 
these fits. The statistical uncertainties on the fit parameters are taken as systematic uncertainties 
("Fit bias"). 

The systematic errors for the parameters that measure interference effects are dominated by 
the uncertainty in the signal model, mainly the tail description of the p resonance. For the other 
parameters, the uncertainty on the fit bias and the S-background contamination are important. 
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Figure 4: Distributions of (clockwise from top left) AE"', rriES) NN output, At/(T(At), m' and 6' for 
samples enhanced in tt^'k~-k^ signal. The dots with error bars give the on-resonance data. 

The solid histogram shows the projection of the fit result. The dark, medium and light shaded 
areas represent respectively the contribution from continuum events, the sum of continuum events 
and the i?-background expectation, and the sum of these and the misreconstructed signal events. 
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Parameter 


Description 


Result 


r 

i_ 


Coefficient of |/— | sin(Am(^Atj 


nio-unii -i-fi no 
— u.iy ± u.ii ± u.Uz 


-'+ 


^~^i~m"ffi r>n1" r»"f 1 ^ . 1 ^ Gin ^ ^ m t ^ 


n n« _|_ n 1 1 _|_ n no 


^0 


Coefficient of 1 fnP 


16 ± 05 ± 05 


U" 


Coefficient of 1 f_P cos(AmwAt) 


0.22 ±0.16 ±0.05 




Coefficient of f _ P 


1.19 ±0.12 ±0.03 




Coefficient of cos(AmdAt) 


0.50 ±0.17 ±0.06 




Coefficient of Imff . f*] cosfAm^At) 


0.3 ±1.4 ±0.3 




Coefficient of Re[/+/*1 cos(AmdAi) 


2.0 ±1.2 ±0.2 




Coefficient of Im[/_|_/1] 


0.16 ± 0.70 ± 0.14 


H — 


Coefficient of Re [/+/*] 


-0.26 ±0.65 ±0.17 




Coefficient of Imf f+ f * 1 sin(AmwAt) 


-5.2 ±1.9 ±0.7 


jRe 


Coefficient of Re[/+/* ] sin(AmdAi) 


-0.3 ±2.0 ±0.5 


7-7+,Im 


Coeflficient of Im[/+/Q] 


0.25 ±0.35 ±0.18 


7-7+, Re 
^+0 


Coefficient of Re[/+/Q] 


-0.34 ±0.39 ±0.15 


r7-+,Im 
•^-0 


Coefficient of Im[/_/Q] 


0.34 ±0.43 ±0.17 


rr+.Re 
^-0 


Coefficient of Re[/_/o*] 


-0.98 ±0.44 ±0.18 



Table 2: Fit results for the U and I coefficients. The errors given are statistical (first) and systematic 
(second). The free normalization parameter is fixed to 1. 



5 FIT RESULTS 

The maximum-likelihood fit results in the vr+vr^vr*^ event yield 1184 ± 58, where the error is 

statistical only. For the U and / coefficients we find the results given together with their statistical 
and systematic errors in Table 2. The corresponding correlation matrix is given in Table 3. We have 
generated a sample of Monte Carlo experiments to determine the probability density distributions 
of the fit parameters. Within the statistical uncertainties of this sample we find symmetric pulls 
and root-mean-squares of order of unity for the U and / coefficients. This allows us to use the 
least-squares method to derive other quantities from these (Section 6). 

The signal is dominated by B° — p^ir^ decays. We observe an excess of /9°7r° events, which 
is in agreement with our previous upper limit of 2.9 x 10~^ at 90% C.L. [22]. The result for the 
p(1450) amplitude is in agreement with the findings in r and e'^e~ decays [11]. For the relative 
strong phase between the p(770) and the /9(1450) amplitudes we find (1991^7)° (statistical error 
only), which is compatible with the result from r and e+e" data. 

Figure 4 shows distributions of Ai?', rriEs, the NN output, At/a{At), where a (At) is the per- 
event error on At, as well as the Dalitz plot variables m' and 9', which are enhanced in signal 
content by cuts on the signal-to-continuum likelihood ratios of the other discriminating variables. 

As a validation of our treatment of the time dependence we allow t^o and Am-d to vary in 
the fit. We find r^o = (1.549 ± 0.080) ps and Am^ = (0.57 ± 0.11) ps~^, while the remaining free 
parameters are consistent with the nominal fit. To validate the SCF modeling, we let the average 
SCF fractions per tagging category free to vary in the fit and find results that are consistent with 
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-0.01 


-0.04 


0.02 


0.01 
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-0.10 


-0.01 
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0.05 


0.13 


0.13 


-0.05 


0.15 


1.00 














-0.03 


0.02 


-0.02 


0.02 
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-0.03 


1.00 












0.09 


-0.05 


0.06 


-0.10 


0.04 
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1.00 








rRo 
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n no 
— U.Uz 


U.UU 


U.Ui 


U.Uo 


U.Uo 
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n no 
— U.Uz 


— U.U4 


U.UO 


—U.Uo 


—U.Ui 


i.UU 






rr+,Im 
^+0 


-0.07 


-0.00 


-0.18 


0.03 


0.03 


-0.01 


-0.11 


0.08 


0.06 


0.02 


-0.10 


0.12 


1.00 




TT+,Re 
^+0 


-0.00 


-0.01 


0.04 


-0.19 


0.01 


0.00 


0.13 


-0.00 


-0.06 


-0.04 


0.13 


0.01 


-0.16 


1.00 


rr+,Im 
^-0 


0.02 


-0.01 


0.13 


0.01 


-0.05 


0.01 


0.03 


-0.01 


-0.04 


0.04 


0.00 


-0.15 


-0.19 


-0.00 1.00 


TT+,Re 
^-0 


-0.05 


-0.05 


-0.05 


0.01 


0.01 


0.04 


-0.06 


0.03 


0.13 


0.01 


-0.07 


0.12 


0.14 


-0.05 -0.01 



Table 3: Correlation matrix for the U and / coefHcients obtained including statistical and systematic 
errors. 



the MC prediction. We also leave the i?-background yields that normalize the S-background classes 
free to vary for those classes that have sufficiently different distributions of the fit variables to be 
distinguished from signal. The results are compatible with the model assumptions. 

6 INTERPRETATION OF THE RESULTS 

The U and I coefficients are related to the quasi-two-body parameters, defined in Ref. [10]^, as 
follows 

^^ = jjf^ ^" = ^' ^^ = ^' ^~ = TF' ^'^ = ^uX+^u+ ' ^^^^ 

and where C = (C+ + C')/2, AC = (C+ - C-)/2, S = {S+ + S-)/2, and AS = (3+ - S-)/2. In 
contrast to our previous analysis [10], the definitions of Eq. (26) explicitly account for the presence 
of interference effects, and are thus exact even for a p with finite width, as long as the U and / 
coefficients are obtained with a Dalitz plot analysis. This treatment leads to a dilution of the result 
and hence to slightly increased statistical uncertainties compared to neglecting the interference 
effects. 

For the CP- violation parameters, we obtain 

Ap^ = -0.088 ± 0.049 ± 0.013 , (27) 
C = 0.34 ± 0.11 ± 0.05 , (28) 
S = -0.10 ±0.14 ±0.04, (29) 



_^For the cosine coefficients we adopt the convention: C+ = [r(B" -> p+7r") -r(B^ p+n~)]/[r{B° p+7r") + 
T{B° p+TT-)] and C" = [V{B° ^ p-7r+) - V{B° p-7r+)l/[r(B° ^ p-7r+) + V{B° p-7r+)]. 
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Figure 5: Confidence level contours for the direct CP asymmetries versus Ap^. The shaded 
areas represent la, 2a and 3a contours, respectively. 



where the first errors given are statistical and the second systematic. For the other parameters in 
the quasi- two-body description of the B^{B^) pir decay-time dependence, we measure 

AC = 0.15 ± 0.11 ± 0.03 , (30) 
AS = 0.22 ±0.15 ±0.03 . (31) 

We find the linear correlation coefficients cc,ac = 0.20 and cs,as = 0.14, while all other correlations 
are smaller. The systematic errors are dominated by the uncertainty on the CP content of the 
-B-related backgrounds. Other contributions are the signal description in the likelihood model 
(including the limit on non-resonant tt^tt^tt^ events), and the fit bias uncertainty. 

One can transform the experimentally convenient, namely uncorrelated, direct CP-violation 
parameters C and Ap^ into the physically more intuitive quantities A'^~, Ap^, defined by 

<-i — P ~ 1 ApTY ± C ± Apj^AC ('lo^ 

■^P- - |k+-|2 + i - l + AC±^p^C ' ^ ' 

^ _ \k — 1 ^ ApTT — C — ApTY^C 

~ |k-+|2±1 ~ 1 - AC - Ap^C ' 

where k"' = {q/p){A~ /A^) and k = {q/p){A^ / A~), so that A'^ {^'^) involves only diagrams 
where the p (vr) meson is emitted by the W boson. We find 

A-j,^ = -0.21 ±0.11 ±0.04 , (33) 
A-p+ = -0.47l°;l^±0.06, (34) 

with a correlation coefficient of 0.59 between A'j^ and A'^^ . The confidence level contours including 
systematic errors are given in Fig. 5. The significance, including systematics, for the observation 
of non-zero direct CP violation is at the 2. Oct level. 
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Figure 6: Confidence level functions for 5^ (left) and a (right). Indicated by the dashed horizontal 

lines are the C.L. values corresponding to la and 2a, respectively. 



The measurement of the resonance interference terms allows us to determine the relative phase 

<5+_ = arg {A+*A-) , (35) 

between the amplitudes of the decays p^n^ and p^vr^. Through the definitions (8)- 

(12), we can derive a constraint on (5_| from the measured U and I coefficients^'^ by performing 

a least-squares minimization with the six complex amplitudes as free parameters. In this fit, one 
complex amplitude can be fixed due to an arbitrary global phase and the normalization condi- 
tion = 1, leaving 10 real-valued unknowns. We obtain the confidence level (C.L.) function 
represented by the dashed line in the left hand plot of Fig. 6. The function includes systematic 
errors. 

This result does not require assumptions beyond the ones outlined in the introduction. The 
constraint can be improved with the use of strong isospin symmetry. The amplitudes A'^ represent 
the sum of tree-level and penguin-type amplitudes, which have different CKM factors: the tree- 
level {T'^) B^ p'^ir'^ transition amplitude^^ is proportional to T/ud^f,, while the corresponding 
penguin-type amplitude {P'^) comes with VqdV*^^, where q = u,c, t. Using the unitarity of the CKM 
matrix one can reorganize the amplitudes and obtains [12] 

A'^ = T^e-'" + P"" , 

^ 7^«e+i« + p^^ (36) 

where the magnitudes of the CKM factors have been absorbed in the T'^, P'^, T'^ and P'^. The 
Eqs. (36) represent 13 unknowns of which two can be fixed due to an arbitrary global phase and the 
normalization condition = 1. Using strong isospin symmetry and neglecting isospin-breaking 
effects, one can identify P^ = —{P'^ + P^)/2, which reduces the number of unknowns to be 
determined by the fit to 9. This set of parameters provides the constraint on 6-\ , given by the 

^"Using directly the relation 6+- = — arctan [((y^l"™ + (7^1^™) / (C/^I^° + !7^I_^°)] to determine 6-f , leads only 

to a weak constraint, because it does not include all the available fit information. 
^^We denote by K the charge conjugate of k, where = 0. 
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solid line on the left plot of Fig. 6. The fit returns a minimum ^ of 7.4, which has a significance 
level of 0.39 for 7 degrees of freedom. We find for the solution that is favored by the fit 



<5+_ = (-6711? ±7)° , 



(37) 



where the first errors are statistical and the second systematic. Only a marginal constraint on 

is obtained for C.L. < 0.05. 

Finally, following the same procedure, we can also derive a constraint on a from the measured 
U and / coefficients. The resulting C.L. function versus a is given in the right hand plot of Fig. 6. 
It includes systematic uncertainties. Ignoring the mirror solution at a + 180°, we find 



where the systematic error is dominated by the uncertainties in the signal model. Only a marginal 
constraint on a is obtained for C.L. < 0.05. 



We have presented the preliminary measurement of CP-violating asymmetries in Br — > 7r''"7r~7r 
decays dominated by the p resonance. The results are obtained from a data sample of 213 millions 
T(4S') BB decays. We extend the previous quasi-two-body approach by using a time-dependent 
Dalitz plot analysis. From the measurement of the coefficients of 16 form factor bilinears we 
determine the three CP-violating and two CP-conserving quasi-two-body parameters, where we 
find a 2.9(7 evidence of direct CP violation. Taking advantage of the interference between the 
p resonances in the Dalitz plot, we derive constraints on the relative strong phase between P" 
decays to p~7r'^ and p^n^ , and on the angle a of the Unitarity Triangle. These measurements are 
consistent with the expectation from the CKM fit [23] . 
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